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Abstract 

We prove analytical results showing that decoherence can be useful for mixing 
time in a continuous-time quantum walk on finite cycles. This complements the nu- 
merical observations by Kendon and Tregenna (Physical Review A 67 (2003), 042315) 
of a similar phenomenon for discrete-time quantum walks. Our analytical treatment of 
continuous-time quantum walks includes a continuous monitoring of all vertices that 
induces the decoherence process. We identify the dynamics of the probability distri- 
bution and observe how mixing times undergo the transition from quantum to classical 
behavior as our decoherence parameter grows from zero to infinity. Our results show 
that, for small rates of decoherence, the mixing time improves linearly with decoher- 
ence, whereas for large rates of decoherence, the mixing time deteriorates linearly 
towards the classical limit. In the middle region of decoherence rates, our numeri- 
cal data confirms the existence of a unique optimal rate for which the mixing time is 
minimized. 



1 Introduction 

The study of quantum walks on graphs has gained considerable interest in quantum com- 
putation due to its potential as an algorithmic technique and as a more natural physical 
model for computation. As in the classical case, there are two important models of quan- 
tum walks, namely, the discrete-time walks 10 |22l El El, and the continuous -time walks 
lfTTT l9ll8l fT0l . Excellent surveys of both models of quantum walks are given in lfT8lfT9ll . In 
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this work, our focus will be on continuous-time quantum walks on graphs and its dynamics 
under decoherence. 

Some promising non-classical dynamics of continuous-time quantum walks were 
shown in ll23l l20l l8l. In ll23ll . Moore and Russell proved that the continuous-time quan- 
tum walk on the n-cube achieves (instantaneous) uniform mixing in time 0(ri), in contrast 
to the fi(nlogn) time needed in the classical random walk. Kempe [ 20 1 showed that the 
hitting time between two diametrically opposite vertices on the n-cube is n°^ l \ as opposed 
to the well-known Vt{2 n ) classical bound (related to the Ehrenfest urn model). In an 
interesting algorithmic application of a continuous-time quantum walk on a specific black- 
box search problem was given. This latter result relied on the exponentially fast hitting 
time of these quantum walks on path-collapsible graphs. 

Further investigations on mixing times for continuous-time quantum walks were given 
in 0[T3|[T1. These works prove non-uniform (average) mixing properties for complete 
multipartite graphs, group-theoretic circulant graphs, and the Cayley graph of the symmet- 
ric group. The latter graph was of considerable interest due to its potential connection to 
the Graph Isomorphism problem, although Gerhardt and Watrous's result in lfT5ll strongly 
discouraged natural approaches based on quantum walks. All of these cited works have 
focused on unitary quantum walks, where we have a closed quantum system without any 
interaction with its environment. 

A more realistic analysis of quantum walks that take into account the effects of deco- 
herence was initiated by Kendon and Tregenna [21]. In that work, Kendon and Tregenna 
made a striking numerical observation that a small amount of decoherence can be useful to 
improve the mixing time of discrete quantum walks on cycles. In this paper, we provide 
an analytical counterpart to Kendon and Tregenna's result for the continuous-time quan- 
tum walk on cycles. Thus showing that the Kendon-Tregenna phenomena is not merely an 
artifact of the discrete-time model, but suggests a fundamental property of decoherence in 
quantum walks. Recent realistic treatment for the hypercube was provided in a recent work 
by Alagic and Russell 0. Developing algorithmic applications that exploit this positive 
effect of decoherence on quantum mixing time provides an interesting challenge for future 
research. 

In this work, we prove that Kendon and Tregenna's observation holds in the continuous- 
time quantum walk model. Our analytical results show that decoherence can improve the 
mixing time in continuous-time quantum walk on cycles. We consider an analytical model 
due to Gurvitz JT31I that incorporates the continuous monitoring of all vertices that induces 
the decoherence process. We identify the dynamics of probability distribution and observe 
how mixing times undergo transition from quantum to classical behavior as decoherence 
parameter grows from to oo. For small rates of decoherence, we observe that mixing 
times improve linearly with decoherence, whereas for large rates, mixing times deteriorate 
linearly towards the classical limit. In the middle region of decoherence rates, we give 
numerical data that confirms the existence of a unique optimal rate for which the mixing 
time is minimal. 
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2 Preliminaries 



Continuous-time quantum walks are well-studied in the physics literature (see, e.g., lfT2ll . 
Chapters 13 and 16), but mainly over constant-dimensional lattices. It was studied recently 
by Farhi, Gutmann, and Childs lfTTl l9l in the algorithmic context. Let G = (V, E) be an 
undirected graph with adjacency matrix Aq. The Laplacian of G is defined as £ = Ac — D, 
where D is a diagonal matrix with Djj is the degree of vertex j 1 . If the time-dependent 
state of the quantum walk is then, by the Schrodinger's equation, we have 



(i) 



The solution of the above equation is \ip(t)) = e~ ltc \vjj(0)) (assuming H = 1). 

We consider the N- vertex cycle graph Cn whose adjacency matrix A is a circulant 
matrix. The eigenvalues of A are Aj = 2 cos(2ir j/N) with corresponding eigenvectors 
\vj), where (k\vj) = H= exp(— 2nijk/N), for j = 0, 1, . . . , N — 1. So, if the initial state 



of the quantum walk is |^(0)) = |0), then \ip(t)) 
terms of the eigenvectors \vj), we get 



e ltL \0). After decomposing |0) in 



\m) 



N-l 



j=0 



-itAj 



Vj). 



(2) 



The scalar term e is an irrelevant phase factor which can be ignored. 



Figure 1: Continuous-time quantum walk on the cycle C 5 . This is a plot of |^ (^)| 2 f° r 
t E [0, 500]. It exhibits a short-term chaotic behavior and a long-term oscillatory behavior. 

If \ip(t)) represents the state of the particle at time t, let Pj(t) = \(j\ip(t))\ 2 be the 
probability that the particle is at vertex j at time t. Let P(t) be the (instantaneous) proba- 
bility distribution of the quantum walk on G. The average probability of vertex j over the 

'We have D — kl, if G is fc-regular. 
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time interval [0,T] is defined as by Pj(T) = ± J* Pj{t) dt. Let P(T) be the (average) 
probability distribution of the quantum walk on G over the time interval [0, T]. 

To define the notion of mixing times of continuous-time quantum walks, we use the 
total variation distance between distributions P and Q that is defined as \\P — Q\\ = 
J2s \P( S ) ~ Q( s )\- F° r £ > 0, the e-mixing time T mix (e) of a continuous-time quan- 
tum walk is the minimum time T so that \\P(T) — Ug\\ < s, where Uq is the uniform 
distribution over G, or 



N-l 



T mix (e) = min < T : 



3=0 



PAT) - v 



< e 



(3) 



Gurvitz's Model To analyze the decoherent continuous-time quantum walk on C N , we 
use an analytical model developed by Gurvitz QUI. In this model, we consider the 
density matrix pit) = \ip(t))(ip(t)\ and study its evolution under a continuous monitoring 
of all vertices of C n . Note that in this case, the probability distribution P(t) of the quantum 
walk is specified by the diagonal elements of pit), that is, Pj(t) = Pj,j{t). 

The time-dependent non-unitary evolution of p(t) in the Gurvitz model is given by (see 
I): 



Pj,k+1 ~~ Pj+l,k — Pj~l,k + Pj,k-1 



r (i - <\jj, ) Pj,k 



Our subsequent analysis will focus on the variable Sj t k defined as 

Sj t k = i J Pj,k 



(4) 



(5) 



The above substitution reduces the system differential equations with complex coefficients 
into the following system with only real coefficients: 



d 1 



Sj,k-i) — r (i — Sj t k) Sj t k- 



(6) 



Throughout the rest of this paper, we will focus on analyzing Equation © for various rates 
of T. One can note that, if T = 0, there is an exact mapping of the quantum walk on a cycle 
onto a classical random walk on a two-dimensional torus. If T ^ 0, there is still an exact 
mapping of the quantum walk on a cycle onto some classical dynamics on a directed toric 
graph. This observation may be useful in estimating quantum speedup in other systems. 



3 Small Decoherence 

We consider the decoherent continuous -time quantum walks when the decoherence rate T 
is small. More specifically, we consider the case when TN 1. First, we rewrite © as 
the perturbed linear operator equation 

±S(jt) = (L + V)S(t), (7) 



4 



Figure 2: The classical recurrence on the 2-dimensional torus derived from a decoherent 
continuous-time quantum walk on the cycle. 

where the linear operators L and U are defined as 

= 4 {^oi,ti^p,v-l + 8 a ^-\8p jV — 5 a)fl 5f3 yU+ i — 8 a)IJ ,+i8p }V ) (8) 

ufci! = -r^ A ,(i-M- (9) 

Here, we consider L as a N 2 x N 2 matrix where L^' v l is the entry of L indexed by the row 
index (//, u) and the column index (a, ft). We view U in a similar manner. The solution of 
© is given by S(t) = e t( ^ + ^S(0), or 

^= E (10) 

fi,v=0 

where < a, /3, fi, v < N — 1. The initial conditions are 

Pa,/3(0) = £ a ,0(O) = £a,0<W (11) 

Perturbation Theory We will use tools from the perturbation theory of linear operators 
(see ifTTlfToll ). To analyze Equation Q, we find the eigenvalues and eigenvectors of L + U. 
Suppose that V is some eigenvector of L with eigenvalue A, that is, hV = XV. Considering 
the perturbed eigenvalue equation 

(L + U)(^ + V) = (A + A) (V + V), (12) 

we drop the second-order terms UV^ and AV^ to obtain the first-order approximation 

VV + hV = XV + XV. (13) 
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By taking the inner product of the above equation with V\ and since L is Hermitian, we 
see that the eigenvalue perturbation term A is defined as 

A = V+UV. (14) 

Let E\ be an eigenspace corresponding to the eigenvalue A and let {Vk : k E 1} be a 
set of eigenvectors of L that spans £\. Let V = J2kei c kVk be a unit vector in E\. Using 
Equation (fl3t . we have 

X)cfcW fc = A^c fc V fc , (15) 
fcei" fee/ 

and after taking the inner product with Vj, we get £) fcg/ c fc^- UV^ = Xcj. If the linear 
combination is uniform, that is Cj = c, for all j, then the eigenvalue perturbation A is simply 
given by 

~ x = J2 v l VVk - (16) 

k&I 

In the case when £\ is one-dimensional or the matrix U is diagonal under all similarity 
actions VjUVk, for j, k <E I, the correction to the eigenvalues is given by the diagonal term 
A = VtUV. Otherwise, we need to solve the system described by det(U\ — XI) = 0. 

To analyze the equation S'(t) = (L + U)S(t), for which the solution is S(t) = 
exp[i(L + U)]5(0), we express 5(0) as a linear combination of the eigenvectors of L + U, 
say {Vj + Vj}. In our case, the evolution of S(t) can be described using the eigenvectors 
of L, since the contribution of the terms Vj are negligible. If 5(0) = J2j c i^v wnere Vj 
are the eigenvectors of L, then 



Spectral Analysis The unperturbed linear operator L has eigenvalues 

/ n(m + ri)\ ( ix(m — n)\ 
\ m , n) = % sin [ K N ' j cos y K N ' ) (18) 

with corresponding eigenvectors 



V M = eX P ( ^ + nZ/ ) ) ' ( 1 9 > 



Thus, for < m, n < N — 1, we have 



N-l 
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To analyze the effects of U, we compute the similarity actions of the eigenvectors on U: 

U(m,rO,(m',„') = (VM)t W ( m '' B ') (21) 

= -^^^-^expt^lim' -m)a + (n -n)b]j (22) 

r 

— — r S m ^ m 8 n i >n + — ^[(m'-m)+(n'-n.)] (mod7V),0 (23) 

where < m, m', n, n' < N — 1. 

The eigenvalues Xr m n ) of L have the following important degeneracies: 

(a) Diagonal (m — n): \{ m>m ) = i sm(2iim/N). 

Each of this eigenvalue has multiplicity 2, by the symmetries of the sine function. 
This degeneracy is absent in our case, since U is diagonal over the corresponding 
eigenvectors. For example, U( m , n ),(N/2-m,N/2-m) =0, for < m < N/2. 

(b) Zero (m + n = (mod N)): A( m>n ) = 0. 

This degeneracy is absent in our case since the corresponding eigenvectors are not 
involved in the linear combination of the initial state S(0). 

(c) Off-diagonal (m ^ n): \ m>n ) = \ n ,m)- 

Since A( TO) „) = % [sin(27rm/A^) + sin(27rm/iV)], each of this eigenvalue has multi- 
plicity at least 4, due to the symmetries of the sine function. In our case, the effective 
degeneracy of these eigenvalues are 2, again by a similar argument. 

By (l23l . the off-diagonal contribution is present if m+n = m'+n' (mod N). Thus, 
A( m ,n) = \m',n') implies that cos(7r(m — n)/N) = ± cos(-7r(m' — n')/N), since 
sin(7r(m + n)/N) = ±sin(7r(m' + n')/N). This implies that m — n = —(m 1 — n') 
or | (m — n) — (m' — n')\ = N, since — (N — 1) < m — n, m! — n' < N — 1. In either 
case, we get m = n ± N/2 or m' = n' ± AT/2. But, upon inspection, we note that U 
is diagonal over these combinations, except for the case when (m', n') = (n, m). 

In what follows, we calculate the eigenvalue perturbation terms A. For simple eigenval- 
ues, these correction terms are given by the diagonal elements 

\m,n) = (V^YVV^ = -r ^" 1 ) , (24) 

by Equation (l23l) . For a degenerate eigenvalue \( m ,n) with multiplicity two, if V = 
c(V r ( m .«) +1/(™- m )), for some constant c, then A( m , n ) = (V^ n ))tUV, and similarly for 
y(n,m)_ Purthej- calculations reveal that the eigenvalue perturbation A( m>n ) is 

A (m ,n) = (y (m>n) ) t U V (m ' n) + (V r(m ' n) ) t U V in ' m) = -T ^' 2 ^ , (25) 
again by Equation (l2"3l . 
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Dynamics We are ready to describe the full solution to Equation ©. First, note that there 
exists a trivial time-independent solution given by S% Jt) = that can be expressed as 
the following linear combination of the eigenvectors of L: 

s°(t) = J2^ + <^+^) y[m,n) - (26) 



N 

(m,n) 

The particular solution will depend on the initial condition 5(0), where S^^O) = S a ,oSpfl- 
Note that we have 



(m,n) 



(m,n) 



Thus, the solution is of the form 



6, 



(m,n) 



2? ™ / n\ 

——{ma + rip) 



(27) 



(28) 



The probability distribution of the continuous-time quantum walk is given by the diagonal 
terms Pj(t) = Sj j(t), that is 



~ A/ + ~m ^ ~ <Wn(modiV 1 n - 



) x S mjn e r n * + (1 - 5 mjn )e 



N 



(m,n) 



x exp 



i sin 



ir(m + n) 
N 



cos 



ir[m — n) 



N 



exp 



2m 

Iv 1 ' 



m + n)j 



We calculate an upper bound on the e-uniform mixing time T mix (e). For this, we define 



iV-l 



it sin(27rm jN)^raj 



N ) 



(29) 



m=0 



where lo n = exp(2iri/N). Note that 



N-l 



N-l 



it\ (m , n) ,,(m+n)j 



m,n=0 

Using these expressions, we have 
1 



(30) 



m=0 



N 



M/(t/2) + 



-tr/jv _ 1 



e -tr/N _ 1 



M 2j (t) 



2-{N mod 2) 



1 + 



iV 



(1 - 2/N) 



(31) 
(32) 



One can note that \Mj{t)\ < 1, and therefore, 

1 



3 

N-l 

E 



^(0 



N 



< e" r — * (JV + e 



-tr/jv 



1). 



(33) 



Since e tT / N < 1, the above equation shows that Ne r n 1 < e. This gives the mixing 
time bound of 

4 Large Decoherence 

We analyze the decoherent continuous -time quantum walks when the decoherence rate T 
is large, that is, when T ^> 1. In our analysis, we will focus on diagonal sums of the matrix 
S(t) from ©. For k = 0, . . . , N — 1, we define the diagonal sum D k as 

AT-l 

, j+k (mod N), 

(35) 

3=0 

where the indices are treated as integers modulo N. We note that 

±D k = -r{l-6 k>0 )D k . (36) 

We refer to the diagonal D as major and the other diagonals as minor. Equation OoT ) 
suggests that the minor diagonal sums decay strongly with characteristic time of order 
1 /r. By the initial conditions, the non-zero elements appear only along the major diagonal. 
From ©, it follows that the system will evolve initially in the following way. The elements 
on the two minor diagonals nearest to the major diagonal will deviate slightly away from 
zero due to nonconformity of classical probability distribution along the major diagonal. 
This process with a rate of order 1/4 will compete with a self-decay with rate of order 
T 3> 1/4, thereby limiting the corresponding off-diagonal elements to small values of the 
order 1/r. A similar argument applies to elements on the other minor diagonals which 
will be kept very small compared to their neighbors that are closer to the major diagonal 
and will be of the order of 1 /T 2 , etc. By retaining only matrix elements that are of order 
of 1 /r, we derive a truncated set of differential equations for the elements along the major 
and the two adjacent minor diagonals: 

= 4 (^i+ij+i _ - rSjj+i, (38) 
/ l 

= 4 ($3,3 ~ Sj-ij-i) - TSjj-!. (39) 
To facilitate our subsequent analysis, we define 
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Then, we observe that 



/ _ ( d j - dj-i) ,i _ (%+i - %) 



a„- = 



d' = v J+ + 3 > - Yd r (41) 



The general solution of the above system of difference equations has the form 

JV-l 

N 



I N ~ 1 

aj = i A k,i ex P (-7fe,o0 + ^fc,2 exp (-7fc.it)} (42) 



fc=0 

1 

d i = { Dk > 1 exp (~Tfc,o*) + Afc,2 exp (-7fc,i*)} ^* 



fc=0 



where c<j = e 2m l N , and the exponents 7 fe)0 and 7 fc) i are the quadratic roots of 



Letting 7 M < 7fc,i. we nave 



1 








sin 2 




2T 








1 




r - 




sin 2 | 




" 2T 





and, for b = 0, 1, we have 



i 



7rk\ ( ink 



(43) 



X (T- X )= 1 -s^(^). (44) 



(45) 



7M = r-^sin^) +0 (i). (46) 
By the initial conditions %(0) = and dj(0) — 0, for j = 0, . . . , N — 1. Thus, 

A M ^ 1, A fc) i ~ " f2 sin2 ^T (47) 



- f sin ^aT ) exp J ' (48) 

These equations show that the amplitudes of the elements along minor diagonals are re- 
duced by an extra factor of T compared to the elements along the major diagonal. Summa- 
rizing, the solution of differential equation at large T has the form 



a i 
3 N 



1 N ~ 1 ( sin 2 \ 

vE ex P[—2r Lt ) uik - (49) 



Based on the above analysis, the full solution for S(t) is given by 



dj if j = k 

Sj,k(t) = { d il 2 if \j -^1 = 1 (50) 
otherwise 
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It can be verified that S(t) is a solution to Equation © modulo terms of order o(l/T). 

The total variation distance between the uniform distribution and the probability distri- 
bution of the decoherent quantum walk on Cjv is given by 



N-l 

E 

i=o 

which simplifies to 

JV-1 



N-l 

E 



JV-1 



-F 



sin 



exp 



2 nk 
N 



k=0 



2T 



t exp 



1 

iV 



E 

j=0 



1 ^ 



JV-l 



^exp 



sin 



2 Trfc 
N 



k=l 



2T 



t cos 



2irkj 
N 



(51) 



(52) 



Lower bound A lower bound on the mixing time for large decoherence rate T can be 
derived as follows. Note that 



N-l 

E 

3=0 



a jit) - 



N 



> 



> 



a (t) - 



N 



k=l 



sin 2 2^ 



2r 



N 



exp 



sin 2 f t 



2r 



(53) 



(54) 



where the first inequality uses the term j = only and the second inequality uses the terms 
k = 1, N — 1. This expression is monotone in t, and is a lower bound on the total variation 
distance. It reaches e at time T[ ower , when 



2T 



- lower 



sin 2 ^ 



In 



2 



2riV 2 



In 



7T Z 



2 



(55) 



for large AT > 1. 



Upper bound An upper bound on the mixing time for large decoherence rate T can be 
derived as follows. Consider the following derivation: 




since | cos(x)| < 1. The last expression is equal to 
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where the last inequality is due to sin(:r) > 2x/-k, whenever < x < n/2 (see Eq. 
4.3.79, Q). Since k > 1, we have k 2 > k. Thus, we have 



E - iv < 2 E exp ( - w J < 2 E ex p I - 

1=0 * k=l V 7 k=l V 



TiV 2 



(60) 



3=0 

The last expression is a geometric series that equals 2/[exp(2£/ (TiV 2 )) — 1]. This expres- 
sion is monotone in t, and it is the upper bound for the total variation distance. It reaches e 
value at time T upper , when 

TN 2 f2 + e\ 

Tupper = — ml— —)■ (61) 

5 Conclusions 

In this work, we studied the average mixing times in a continuous-time quantum walk on 
the iV-vertex cycle Cn under decoherence. For this, we used an analytical model developed 
by S. Gurvitz |[T3l . We found two distinct dynamics of the quantum walk based on the rates 
of the decoherence parameter. For small decoherence rates, where TN <C 1, the mixing 
time is bounded as 

1 . fN s 



T m ix *C j_, In 



1 + 



N-2 



(62) 



This bound shows that T mix is inversely proportional to the decoherence rate T. For large 
decoherence rates r> 1, the mixing times are bounded as 

These bounds are show that T mix is linearly proportional to the decoherence rate T, but is 
quadratically dependent on N. Note that the dependences on N of the mixing times exhibit 
the expected quantum to classical transition. 

These analytical results already point to the existence of an optimal decoherence rate 
for which the mixing time is minimum. Our additional numerical experiments (see Figure 
©) for T ~ 1 confirmed that there is a unique optimal decoherence rate for which the 
mixing time is minimum. This provides a continuous-time analogue of the Kendon and 
Tregenna results in EH . 
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